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The orbits of open clusters in the Galaxy 
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ABSTRACT 

We present and analyze kinematics and orbits for a sample of 488 open clusters in the 
Galaxy. The velocity ellipsoid for our present sample is derived as {crjj, ay ■, crvy)=(28.7, 
15.8, 11.0) km s^^ which represents a young thin disc population. We also confirm 
that the velocity dispersions increase with the age of cluster subsample. The orbits of 
open clusters are calculated with three Galactic gravitational potential models. The 
errors of orbital parameters are also calculated considering the intrinsic variation of 
the orbital parameters and the effects of observational uncertainties. The observational 
uncertainties dominate the errors of derived orbital parameters. The vertical motions 
of clusters calculated using different Galactic disc models are rather different. The 
observed radial metallicity gradient of clusters is derived with a slope of 6 = —0.070 ± 
0.011 dex kpc^^. The radial metallicity gradient of clusters based on their apogalactic 
distances is also derived with a slope of 6 = — 0.082±0.014 dex kpc^^. The distribution 
of derived orbital eccentricities for open clusters is very similar to the one derived for 
the field population of dwarfs and giants in the thin disc. 
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1 INTRODUCTION 

Open clusters (OCs) have long been used as important tools 
in the study of the Galactic disc. The young clusters have 
been used to determine spiral arm structure, and to map the 
rotation curve of the Galaxy. The old clusters are excellent 
tracers of the structure, kinematics, and chemistry of the 
Galactic disc (Friel 199f|) . In principle, basic parameters like 
distance, age, and metallicity can be determined for a cluster 
more accurately than for a field star. Therefore, OCs are 
better tracers of large scale properties of the Ga lactic disc 
population than field stars (|Piskunov et a l.' 2006). 

About 20 years ago, the cat alogue of open cluster data 
compiled by Lynga l|Lyngall 19871 ) was used to derive the ra- 
dial grad ients and other structural properties of th e Galac- 
tic disc (|Lvngal [l98^ : iJanes. Tillev fc Lvngal [T98 g ) . In re- 
cent years, the available data on open clusters have in- 
creased very quickly and the basic p arameters of these clus- 
ters have also considerably improved. iDias et al. I (|2002l . here- 
after DAML) compiled the new catalogue of optically visible 
open clusters and candida tes including 1. 5 times more clus- 
ters than the catalogue of lLvngal (| 19871 ). iKharchenko" et al.l 
l|2005l . hereafter K05) presented a catalogue of astrophys- 
ical data for 520 OCs, with data derived fr om the All- 
Sky Compiled Catalogue of 2.5 Million Stars (|Kharchenkd 
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l200ll . hereafter ASCC-2.5). Using these data, the prop- 
erties of the Galactic disc such as the scale height, the 
scale length, and the metallicity distribution of the disc , 
were derived (|Chen. Hou fc Wang|2003l : lBonatto et al ]|2006l : 



iPiskunov et al.ll2006l ). 



iLynga fc Palouj (|l987l ) used 106 OCs with available 
data on positions, distances, and radial velocities to an- 
alyze the local kinematics. They found that the disper- 
sion of the radial velocity increases with age: for the 
old clusters it is about twice that of young clusters. 
iBarkhatova. Kutuzov fc Osipkovl l| 19871 . hereafter BKO) cal- 
culated the Galactic orbits for 69 OCs based on care- 
fully selected data o n distances, absolute proper mot ions, 
and radial velocities. iKeenan. Innanen fc Housa (|l973l ) cal- 
culated the orbits of NGC 188 and M67. In addition, 
lAllen fc MartosI l|l988l ) calculate d the orb its for NGC 188, 
M67, and NGC 2420, [Carraro fc Chiosil (jl994 . hereafter 
C94) expanded th e sample with 5 classic clusters, and 
iFinlav et all (Il995l) calculated the orbits for a total of 7 
old clusters. ISoubiran. Odenkirchen fc Le CampionI 
calculated the orbit for NGC 2355, de Oliveira et al] 



2001 



200: 



calculated the o rbits for NGC 1912 and NGC 1907, and 
IWu et al.l (|2002l ) calculated the orbit for M48. 



In recent years, iDias fc Lepin3 |200i) determined the 
rotation velocity of the spiral pattern of the Galaxy by 
studying the birthplaces of OCs in the Galactic disc as 
a function of their ages. The birthplaces of these clusters 
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were determined by assuming that their orbits were 
circular. Using 148 OCs witliin the projected distan c e onto 
the Galactic plane dxy s; 0.85 kpc. IPiskunov et all (|2006l ) 
derived the Solar motion and the velocity ellipsoid of OCs. 
They also calculated the Galactic orbits of these clusters and 
presented the mean parameters of their orbi ts. Based on the 
spatia l and velocity distributions of OCs, IPiskunov et al.l 
(|2006l l also identified the existence of four open cluster 
complexes of diff erent ages (|Kharchenko fc Piskunovll2006l : 
iRoser et al.ll2007h . which verifed the nature of cluster- 
ing in OCs in the Solar neighbourhood i dentified in 
previ ous studie s (Efremov 1978; Eigens on Si Yatsvkl 



19881; 



Barkhatova. Osipkov &£ Kutuzov 



1989l V 



Lepine. Dias fc Mishurov, (,2008 ) used 374 OCs taken from 



the 2.7 version of DAML catalogue to measure the epicycle 
frequency k in the Galactic disc. They also calculated the 
orbits of these OCs and discussed the distribution of initial 
velocities of these clusters. 

The orbital motions of OCs are important not only 
for our unders tanding of the dynamical evolution of OC s 
in the Galaxy (|Friellll999l : iBergond. Leon fc Guibertll200lf ). 
but also for investigating their effect on the time-evolution 
of the abundance gradient in the Galactic disc (C94). The 
main aim of this paper is to calculate the orbits of OCs with 
the improved data for an enlarged sample of 488 OCs and 
to discuss the kinematical properties of these clusters. 

We present the collection of data in Section 2. In Section 
3, we analyze the statistical properties of the sample, and 
especially their age-velocity dispersion relations. We present 
the orbital parameters and their associated uncertainties cal- 
culated in a given Galactic potential model in Section 4, 
followed by Section 5, which presents the differences of or- 
bital parameters due to different Galactic potential models. 
In Section 6, we compare our results with the orbital pa- 
rameters calculated for field stars and globular clusters, and 
discuss the effect of the orbital motions of OCs for the radial 
abundance gradient derived for these clusters, we also com- 
pare the orbital parameters with those derived by previous 
studies. Our conclusions are given in Section 7. 



proper motions and 447 have radial velocities. 869 clusters 
with distances and proper motions data are taken from the 
DAML catalogue as the initial sample. 



2.2 Distances 

Most of distances and ages listed in DAML are taken 

from WEBDA, and are updated with new data from 

I ' ( I 1 

literature (IDias et a l. 2002). The distances estimated by 

iBaumgardt. Dettbarn fc Widen (2000) based on Hipparcos 
parallaxes ||ESA||1997| ) of member stars are also adopted by 
DAML. A comparison of the Hipparcos pa rallaxes with pho- 
tome tric distances shows good agreement (jBaumgardt et ahl 
[2000^. The distances listed in DAML are taken from different 
authors using different observational techniques and reduc- 
tion methods, and no errors for this parameter are listed in 
DAML, the non- uniformity and uncertainties in the adopted 
dist ance data should be estimated. 

IPaunzen fc Netopill (2006) studied the accuracy of avail- 
able parameters such as the age, reddening, and distance for 
OCs by using the independently derived values published in 
the literature. They used a sample of 395 clusters in their 
statistical analysis. They found that, for about 80% clus- 
ters in their sample, the error of distance is less than 20%. 
They compared their results with the data of DAML cata- 
logue, and pointed out that the distances listed in DAML 
catalogue have the same error distributions. In our present 
study, 20% relative errors in distances are adopted for OCs 
in o ur sample. 

Ivan LeeuwenI (|2008l ) re-reduced the raw data of Hippar- 
cos mission, and the new reduction provides an improve- 
ment by a factor 2.2 comp ared to the catalogue published 
in 1997 (Ivan Leeuwen|[2007h . We calculate the difference be- 
tween the distances derived bv lvan LeeuwenI l|2008h using the 
re-reduced Hipparcos data and that listed in DAML cata- 
logue for a sample of 17 OCs in common. We get a difference 
of 6.2 ± 1.6% between these two data sets, which is less than 
our adopted 20% error for distance parameters and is con- 
sidered in the following orbit calculation. 



2 THE SAMPLE 

2.1 The DAML catalogue 

We chose the DAML catalogue as the main source of the 
fundamental parameters for the OCs. This catalog ue uses 
the W EBDA databascQ and previous catalogues of iLvngH 
|l983) as a starting point. Kinematical, and meta llicity data 
of th e new objects when available are inserted l|Dias et al.l 
120021 '). They also made use of the Simbad databas43 and of 
the literature to find data on the clusters or on individual 
stars of the clusters, to obtain the averaged values of radial 
velocities and proper motions. This catalogue is regularly 
updated, and the latest version is available from its web- 
sitejfl. The present 2.9 version (13/Apr/2008) of DAML cat- 
alogue contains 1776 objects, of which 936 have published 
distances, ages, and reddening values, 890 have published 



^ http://obswww.unige.ch/webda 
^ http:/ /simbad. u-strasbg.fr/simbad/ 
http: / /www. astro. iag.usp.br/" wilton/. 



2.3 Absolute proper motions 

The absolute proper motion data of OCs listed in DAML 
are adopted in the present study. The proper motion data 
are also compiled from different authors, but they are all 
based on the Hipparcos system. In our final sample, proper 
motio ns of 15% c lusters are derived from the Hipparcos cata- 
logue l|ESAlll997l ). 47% proper mot ions are derived from the 
Tycho2 catalogue (|H0g et al ]|2000l ). 36% proper motions ar e 
derived from the ASCC-2.5 catalogue (|Kharchenkd 120011 ) . 
and only 2% proper motions are derived from the UCAC2 
catalogue (Zacharias et al. 2004]). 

The Hipparcos catalogue is considered as the re- 
alization of the ICRS at optical wavelengths. The 
systematic error of proper motions in the Hippar- 
cos catalogue wit h respect t o ICRS is estimated t o 
be 0.25 mas vr-^ JeSAI Il997l: iKovalevskv et al] Il997l ). 
iPlatais. Kozhurina-Platais fc van Leeuwen l|l998h searched 
the Hipparcos catalogue and found 9 n ew OCs and de- 
rived t he proper motions for these clusters. IBaumgardt et al] 
(|2000l ) determined mean proper motions of 205 OCs from 
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their member stars found in the Hipparcos catalogue. In our 
final sample, the proper motions derive d from Hipparcos cat - 
alogue are al l taken from th e results of lPlatais et al.l (|l998l ) 
and .Baumgardt et al.1 (|200d ). 



The Tycho2 catalogue l|H0g et al.ll2000l ) presents very 
precise proper motions with random errors between 1 and 
3 mas yr~^ in the Hipparcos system. There are no signif- 
icant systematic differe nces between the proper motion s 
of these two catalogues (lUrban. Wvcoff Sz Makarovll200Cf ). 
iDias. Leoine fc Alessil (|200ll . l2002l ) determined the mean ab- 
solute proper motions of 206 OC s from the data in Ty- 
cho2 catalogue. lAlessi et af] l|2003l ') found 11 new OCs candi- 
dates in Tycho2 catalogue a nd determined the mean pr oper 
motions for these clusters. iLoktin fc Beshenovl (|2003h de- 
termined the mean proper motions for 167 OCs based on 
the kinematic and photome t ric da ta in Tycho2 catalogue. 
iDias. Lepine fc Al essi (2001, 2002) compared their results 
with those derived by Baumgardt ct al. ( 2 000l ) based on Hip- 
parcos catalogue, and found that the m ean difference in the 
prope r motions is less than 1 mas vr~^. lLoktin fc Beshenovl 
(120031) also compared their result s with those derived by 
iDias. Lepine fc Alessil (|200ll . l2002h and found that the dif- 
ference is 4 ± 5%. In our final sample, the mean proper mo- 
tions derived from Tycho2 catalogue are taken from above 
mentioned studies. 

The ASCC-2.5 catalogue is based on large, mod- 
ern, high- precision catalogues of the Hipparcos- Tycho fam- 
ily, including the Tycho2 catalogue, and provides the 
most complete all-sky catalogue of stars having uni- 
form high precision astrom etric and photom etric data 
down to ~ 14 mag (|Kharchenkol 120011 ). In our 
final sample, the me an proper motions der i ved from 
ASCC-2. 5 catalogue by iKharchenko et al.1 (|2003l . |2005| ) are 
adopted. IKharchenko et al.1 (|2003|. |2005|) compa r ed the ir re- 
sults with those deri v ed by iBaumgardt et al.l (|2000t ) and 
iDias. Lepine fc Alessil (|200ll . 120021 ) and found that their re- 
sults agree quite well with previous studies. 

The UCAC2 catalogue presents proper motions in the 
Hipparcos system with nominal errors of 1 to 3 mas yr~^ 
for stars up to 1/ ~ 12 mag and about 4 to 7 mas yr~^ 
for fainter stars up to 1^ ~ 16 mag. The systematic er- 
rors of the prope r motions in UCAC2 ar e in the ra, n ge 0.5 
to 1.0 mas yr"^ (|Zacharias et al.ll2004h . iDias et al.1 ()2006l ) 
determined the mean proper motions for 428 OCs from 
the UCAC2 catalogue. They compared their res ults with 
those derived from Hipparcos catalogue (.Baumgardt et al.l 
20001) , Tycho2 catalogue (iDias. Lepine fc Alessil l200ll . 
20021: iLoktin fc Beshen"wl |2003| ) and ASCC-2.5 catalogue 
( Kharchenko et al.ll2003l ). and found that there is no statisti- 
cal distinction between the compared mean proper motions 
of OCs in these catalogues. 

We also calc ulate the difference o f the mean proper mo- 
tions derived by Ivan LeeuwenI (|2008l ) using the re-reduced 
Hipparcos data and that adopted in our final sample taken 
from the DAML catalogue for the 17 OCs in common. We 
get a difference of 0.35 ± 0.22 mas yr~^ in /i^cos^ and a 
difference of 0.46 ±0.43 in ns- The differences are within the 
range of the errors for the mean proper motions listed in the 
DAML catalogue. 



2.4 Radial velocities 

There are 431 OCs whose distances, proper motions, and ra- 
dial velocities are available from the DAML catalogue. Most 
ra dial velocities a re taken from the catalogues compiled 
by lKharchenko e t al. (2005, 2007). Kharchenko ct al. (200i) 
cross-identified the ASCC-2.5 cat alogue with the General 
Catalogue of Radial Velocities of iBarbier-Brossat fc FigonI 
(|2000h and derived radial velocities for 290 OCs based on 
their membership determination. They found that the mean 
difference between their results and those in the published 
literature for common clusters is 0.36 ± 0.88 km s"^ The 
363 ra dial velocities of OCs compiled bv IKharchenko et al.l 
(|2007l ) are derived from the 2nd version of the Catalogue of 
Radial Velocities of Galactic stars with high precision As- 
trometric Data (CRVAD-2). The CRVAD-2 is the resuh of 
updating and expanding the list of stars with known radial 
velocities and high precision astrometric and photometric 
data taken from the ASCC-2.5 catalogue. The mean differ- 
ence of radial velocities for 177 clusters in common between 
the CRVAD-2 and the literat ure is 0.6 5 ± 0.72 km s"\ 

iFrinchabov fc Maiewskil (|2008h presented high- 
precision radial velocities for 71 OCs obtained with 
multi-object spectrographs. For 25 clusters in their sample, 
the radial velocities are newly obtained and a re not included 
in DAML. iMermiUiod. Mavor fc Udrvl (|2008l ) derived mean 
radial velocities for 166 OCs based on observations with 
the CORAVEL spectrovelocimeters. The radial velocities of 
64 OCs derived by JMormilliod ct al. (2008) are not listed in 
DAML. We supplement t he ne w ra dial velocities derived h 
IFrinchabov fc Maiewskil (|2008l ) and lMermiUiod et all (|200; ' 
into our final sample. We also update the known radial 
velocities in DAML with those derived bv IMermilliod et al.l 
( 2008. ) due to their higher precision. 

In our final sample, 488 OCs with distances, proper mo- 
tions, and radial velocities are included. 



2.5 Completeness of present sample 

Our present sample is a subsample mainly taken from the 
DAML catalogue. We plot the volume density of 1082 OCs 
(open circles) with distance data available in DAML as a 
function of heliocentric distance in Fig. [1] In order to es- 
timate the completeness of the OCs in our present sample, 
we also plot the volume density of 488 OCs (filled circles) 
in our present sample in Fig. [T] For clusters with heliocen- 
tric distance > 2.0 kpc, the volume density distributions 
for these two samples of OCs are completely consistent. For 
clusters with heliocentric distance d© < 2.0 kpc, the volume 
densities of OCs in our present sample are less than those 
for OCs in DAML, but the distributions for these two sam- 
ples are very similar. Fig. [1] indicates that our present OCs 
sample is a representative subsample of currently observed 
OCs in the Galaxy. 

Using a sample o f 654 OCs with distance data available, 
iBonatto et al.l (|2006l ) simulated the effects of completeness 
in their OCs sample. They found that a total number of 
~ 730 OCs with heliocentric distance cIq ^ 1.3 kpc should 
be observed. Within the same distance range, there are 498 
OCs in DAML catalogue and 271 OCs in our present sam- 
ple, the completeness can be estimated as 68% and 37% for 
DAML catalogue and for our present sample, respectively. 
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A possible source for the remainder unobserved OCs comes 
from the very young clusters which are still embedded within 
giant molecular clouds, they are heavily ob scured and are 
very difficult to identify l|Lada fc Ladal 120031 ). 



2.6 Initial conditions 

Table [1] lists the adopted positions {a,S), heliocentric 



distances dQ, radial velocities 



and absolute proper 
motions {fj,a cos 5, jj,s) of 488 OCs in our final sam- 
ple. The adopted 20% errors for distances, the er- 
rors of absolute pro per motions and of radial v eloci- 
ties listed in DAML, iFrinchabov fc Maiewskil (|2008h . and 
iMermilliod et al ] (|2 008 l) are also presented in Table [l] In Ta- 
ble [TJ we also list the ages and metallicities [Fe/H] for each 
cluster if available in DAML. Ages are currently available 
for 445 OCs, 109 of which have [Fe/H] values also available. 

The initial conditions for orbit calculation are the 
presently observed positions and velocities of OCs with re- 
spect to the galactocentric reference frame. Adopting the 
solar motion (U, V W) r^ = (10.0, 5.2, 7.2) km s~^ from 
iDehnen fc BinnevI (119981 ). the local standard of rest (LSR) 
velocities of OCs are determined from the data listed in Ta- 
ble [T] The LSR velocities are then corrected to the Galactic 
standard of rest (GSR) by a dopting the galactocentric dis- 
tance of Sun Rq = 8.0 kpc (lReidlll993l) and a rotation ve - 
locity of the LSR of 220 km s~^ (|Kerr fc Lvnden-Belill 19861 ). 
All space coordinates x, y, z, and velocity components U,V, 
W refer to a galactocentric righthanded cartesian coordinate 
system with the x direction directed towards the Galactic 
anticenter and the z direction directed to wards the Galactic 
north pole ("Odenkirch en fc Broschll 19921 ). The U,V,W ve- 
locity com ponents and their error s are c alculated with the 
method of IJohnson fc Soderbloinl (|l987l ). Table [5] lists the 
initial conditions used to calculate the orbital solution for 
each cluster in our sample. The errors in velocity compo- 
nents include the errors in absolute proper motions, radial 
velocities, and distances of OCs listed in Table [T] 



3 PARAMETER DISTRIBUTIONS AND 
KINEMATICS OF THE SAMPLE 

3.1 Parameter distributions of the present sample 

Fig. [2] shows the distributions of ages, metallicities [Fe/H], 
and observed galactocentric distances Rqq of OCs in our 
present sample. Due to the large difference among the ages of 
OCs, the distributions of age are plotted in panels a and b of 
Fig.[2l We can see from panels a and b of Fig.[2]that, about 
2/3 OCs have ages less than 100 Myr. The oldest cluster has 
an age of 9.0 Gyr. No clusters in our sample have been found 
to be in the age interval 0.9 - 1.0 Gyr. The distribution of 
metallicities [Fe/H] of OCs is plotted in panel c of Fig. [5] Ex- 
cept for two clusters with [Fe/H] < —0.5, the distribution of 
[Fe/H] can be fitted by two Gaussian functions. For clusters 
with -0.5 < [Fe/H] < -0.2, the [Fe/H] data can be fitted 
by a Gaussian function with mean Mjpgy'jj] = —0.31 and 
dispersion ''"[pg^jj] = 0.07; for clusters with [Fe/H] > —0.2, 

■i] = 0-0 

0.13. It should be noted that the 



the best-fitting Gaussian function has mean fi 



and dispersion a 



[Fe/H] 



sample of clusters with [Fe/H] data is very incomplete, the 
distribution of [Fe/H] in panel c of Fig.[2]may not be the true 
distribution of the metallicities for the complete OCs sample 
in the Galaxy. From panel d of Fig. [2] it can be seen that 
most OCs distribute near the Sun. The minimum and max- 
imum galactocentric distances of OCs in our present sample 
are 4.6 and 22.6 kpc respectively. 

3.2 Age-velocity dispersion relations 

In recent years, many observational efforts have been de- 
voted to c onstrain the age-velocity dispersion relation of the 
thin disc. INordstrom et all l|2004l . hereafter N04) presented 
new determinations of metallicity, age, kinematics, and 
Galactic orbits for a complete, magnitude-limited, and kine- 
matically unbiased sample of ~ 14000 F and G dwarf stars 
near the Solar neighborhood. The Hipparcos/Tycho-2 par- 
allaxes and proper motions, together with some 63000 new, 
accurate radial velocity observations supplemented by a few 
earlier radial velocities, were used to compute the space ve- 
locity components and their dispersions. Ages and their er- 
rors were computed from a set of t heoretical isochrones by 
a sop histicated Bayesian technique (|j0rgensen fc LindegrenI 
l2005h . N04 found that the age- velocity dispersion relations 
of each space velocity component can be fitted by continu- 
ous smooth power laws: a cc age*, which also evidences the 
con tinuous heating of the disc in all dir ect ions. 

iHolmberg. Nordstrom fc Andersei] (|2007l ) redeter- 
mined the basic calibrations used to infer astrophysical 
parameters for the N04 stars from uvby photometry. Using 
the improved astrophysical parameters, they recomputed 
the ages and age error estimates for the N04 sample. Based 
on their revised data set, and with substant ially higher time 
resolution than that in the original N04, IHolmberg et al.l 
(2007) confirmed the conclusion of N04 that the dynamical 
hea ting of the thin disc continu es throughout its life. 

ISeabroke fc Gilmord (|2007l ) revisited the Galactic thin 
disc age-velocity dispersion relation based on the N04 sam- 
ple, their new result is that a power law is not required by 
the data of N04, and disc heating models that saturate after 
~ 4. 5 Gyr are equa ll y con sistent with the observations. 

ISoubiran et al.l 1 20081 . hereafter S08) presented the pa- 
rameters of 891 stars, mostly local and distant clump giants, 
including distances, absolute magnitudes, spatial velocities, 
galactic orbits, and ages. Using their distant sample of clump 
giants, and rejecting stars having a probability higher than 
80% to belong to the thick disc, the Hercules stream, and 
the halo, S08 found that the velocity dispersions in V and 
W saturate at 4 Gyr and the dispersion in U increases 
smoothly with time. 

Usin g radial velocitie s of 67 clusters within 2 kpc from 
the Sun, iLvnga fc Palou^ (|l987h found that the dispersion 
of the radial velocity increases with age. Based on proper 
motions and distances of 148 clusters within the projected 
dista nce onto the G a lactic plane d^y ^ 0.85 kpc from the 
Sun, IPiskunov et al.l (|2006l ) derived the tangential velocity 
dispersions for clusters with different ages and also found 
that the dispersions increase with age. 

The OCs in our sample with errors in the spatial ve- 
locities less than 20 km s~^ are used to derive the velocity 
dispersion in each velocity component. Those clusters are 
divided into three age groups: < age ^ 500 Myr, 500 < 
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Table 1. The observed data of 488 OCs in our sample. The full version of this table is available in the online version of this article. 



name 


a(2000.0) 


(5(2000.0) 






fia cos 5 


Mi 


Age 


[Fe/H] 




h 


m 


s 





/ // 




(kpc) 


(km s~^) 


(mas yr"-"-) 


(mas yr"-"-) 


(Myr) 




Berkeley 59 


00 : 


02 


: 14 


+67 : 


25 : 


00 


1.000 + 0.200 


-12.5 + 7.1 


-2.11 + 0.81 


-1.20 + 0.75 


6.3 




Blanco 1 


00 : 


04 


; 07 


-29 : 


50 : 


00 


0.269 ± 0.054 


4.1 + 1.4 


20.17 + 0.51 


3.00 + 0.51 


62.5 


0.04 


Alessi 20 


00 : 


09 


: 23 


+58 : 


39 : 


57 


0.450 ± 0.090 


-11.5 + 0.0 


8.73 + 0.53 


-3.11 + 0.53 


166.0 




ASCC 1 


00 : 


09 


: 36 


+62 : 


40 : 


48 


4.000 ± 0.800 


-69.7 + 4.7 


-2.07 + 0.72 


0.46 + 0.57 


177.8 




Mayer 1 


00 : 


21 


; 54 


+61 : 


45 : 


00 


1.429 ± 0.286 


-20.9 + 2.0 


-4.46 + 1.13 


-6.66 + 0.94 






NGC 129 


00 : 


30 


; 00 


+60 : 


13 : 


06 


1.625 ± 0.325 


-39.4 + 0.5 


-1.06 + 0.94 


1.60 + 0.94 


76.9 




ASCC 3 


00 : 


31 


; 09 


+55 : 


16 : 


48 


1.700 ± 0.340 


-37.0 + 0.0 


-1.92 + 0.61 


-1.25 + 0.59 


79.4 




NGC 225 


00 : 


43 


: 39 


+61 : 


46 : 


30 


0.657 + 0.131 


-28.0 + 0.0 


-4.95 + 0.76 


-0.50 + 0.76 


130.0 




NGC 188 


00 : 


47 


: 28 


+85 : 


15 : 


18 


2.047 ± 0.409 


-45.0 + 10.0 


-1.48 + 1.25 


-0.56 + 1.24 


4285.5 


-0.01 


IC 1590 


00 : 


52 


: 49 


+56 : 


37 : 


42 


2.940 + 0.588 


-32.5 + 6.4 


-1.36 + 0.23 


-1.34 + 0.83 


3.5 





Table 2. The present positions and velocities of 488 OCs in our sample. The full version of this table is available in the online version 
of this article. 



name 


X 




y 


z 


U 


V 




W 










(kpc) 






(km s~ 






Berkeley 59 


-8.471 


± 


0.094 


0.878 + 0.176 


0.087 + 0.017 


25.3 + 5.1 


219.7 + 


6.6 


2.4 + 0.7 


Blanco 1 


-7.952 


± 


0.010 


0.013 + 0.003 


-0.264 + 0.053 


-13.4 + 4.9 


217.1 + 


1.8 


-1.7 + 0.7 


Alessi 20 


-8.207 


± 


0.041 


0.398 ± 0.080 


-0.030 + 0.006 


0.2 + 3.2 


206.5 ± 


1.8 


-1.6 + 0.2 


ASCC 1 


-9.887 


± 


0.377 


3.527 + 0.705 


0.014 + 0.003 


75.8 + 13.9 


181.3 ± 


8.4 


21.9 + 0.3 


Mayer 1 


-8.702 


± 


0.141 


1.244 + 0.249 


-0.023 + 0.005 


51.2 + 9.1 


223.7 + 


5.3 


-33.8 + 0.5 


NGC 129 


-8.818 


± 


0.164 


1.402 + 0.280 


-0.072 + 0.014 


35.7 + 6.4 


195.3 ± 


3.9 


21.9 + 0.7 


ASCC 3 


-8.843 


± 


0.169 


1.459 + 0.292 


-0.221 + 0.044 


43.0 + 5.2 


200.6 ± 


3.0 


3.3 + 0.0 


NGC 225 


-8.348 


± 


0.069 


0.557 + 0.111 


-0.012 + 0.002 


38.0 + 3.3 


209.6 ± 


2.1 


6.6 + 0.4 


NGC 188 


-9.027 


± 


0.205 


1.590 + 0.318 


0.780 ± 0.156 


43.6 + 11.8 


199.7 + 


1.0 


-14.7 + 0.8 


IC 1590 


-9.597 


± 


0.319 


2.448 ± 0.490 


-0.320 + 0.064 


44.6 + 5.6 


206.7 + 


6.0 


-7.9 + 0.1 



age < 1000 Myr, and 1000 < age 2000 Myr. For each age 
group, the velocity dispersion in each velocity component U, 
V, and W is calculated. The number A'^ of clusters in each 
age group is also listed in Table |3l Table [3] indicates that the 
velocity dispersions in the U and W components in the age 
group of 500 < age ^ 1000 Myr are only marginally larger 
than those in the age group with age < 500 Myr, but the 
velocity dispersions in all of the three velocity components 
in the age group vsrith age > 1 Gyr are larger than those for 
the clusters with age < 1 Gyr. The velocity dispersions in 
our present OCs sample indicate the continuous dynamical 
heating of the thin disc. 

Using OCs with errors in the spatial velocities less than 
20 km s~^ in our present sample, the velocity dispersions 
are derived as (cry, av, <yw) =(28.7, 15.8, 11.0) km s 
The derived velocity dispersions are smaller than those for 
the thin disc clump giants (ajj, av, aw) ~ (41.5, 26.4, 22.1) 
km s ~^ derived by S08. For their subsample of clump giants 
with age of 1.5 Gyr, S08 derived the velocity dispersions as 
((7(7, av, aw) ~ (36.2, 19.9, 18.7) km s which are close 
to our results for the OCs with age > 1 Gyr. So, the main 
reason for the difference between the velocity dispersions for 
OCs in our present sample and those for the giants of S08 
is that in our present OCs sample, most of them are young 
clusters with age less than 1 Gyr. The giants sample of S08 
represents an old thin disc, most of giants in their sample 



are older than 1 Gyr. Just as we have pointed out in the 
previous section, the thin disc of the Galaxy is continuously 
heated, so the velocity dispersions for a young population 
such as the OCs in our present sample should be less than 
that corresponding to an older population such as the giants 
sample of SOS. 



Using OCs within the projected dis tance onto the 
Galac tic plane d^y ^ 0.85 kpc from the Sun. lPiskunov et al.l 
(|2006t ) derived the velocity dispersions as [au, av, 
aw) =(13.86, 8.75, 5.05) km s whi ch are smaller 
than those derived from our present sample. IPiskunov et al.l 
(2006) used a volume-limited sample of relatively young 
clusters; in contrast, our present sample includes many older 
clusters and also objects located further away. Therefore, 
our larger velocity dispersions reflect the effects of dynami- 
cal heating of the thin disc. Using the clusters in our present 
sample within the projected distance onto the Galactic plane 
dxy ^ 0.85 kpc, we derive the velocity dispersions as {au, 
av, aw) =(16.8, 9.6, 5.3 ) km s which are close to the 
results of IPiskunov et all l|2006l ). 
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Table 4. The orbital parameters and their errors of 488 OCs in our sample calculated with AS91 model. The full version of this table 
is available in the online version of this article. 



name 


Ra 


Rp 


e 


•2 m ax 






J. 




(kpc) 


(kpc) 




(kpc) 


Myr 


Myr 


(kpc km s^-*-) 


Berkeley 59 


8.6 ±0.4 


8.5 ±0.3 


0.01 ±0.03 


0.09 ±0.02 


239.3 ±8.4 


35.1 ± 1.6 


-1883.3 ± 65.2 


Blanco 1 


8.2 ±0.1 


7.5 ±0.2 


0.05 ±0.02 


0.26 ±0.05 


221.1 ± 1.3 


36.4 ± 1.3 


-1726.2 ± 16.8 


Alessi 20 


8.3 ±0.1 


7.2 ±0.1 


0.07 ±0.02 


0.03 ±0.01 


217.0 ±0.6 


30.8 ±0.3 


-1694.8 ±8.5 


ASCC 1 


10.6 ±0.8 


8.5 ±1.2 


0.11 ±0.05 


0.33 ±0.26 


268.5 ±27.3 


46.5 ± 11.9 


-2059.9 ±209.0 


Mayer 1 


10.0 ±0.9 


8.5 ±0.3 


0.08 ± 0.03 


0.50 ±0.24 


261.8 ± 16.5 


50.0 ± 10.1 


-2010.3 ±94.4 


NGC 129 


8.9 ±0.2 


7.4 ±0.4 


0.10 ±0.02 


0.28 ±0.12 


229.3 ±7.8 


38.5 ±4.5 


-1772.2 ± 64.7 


ASCC 3 


9.0 ±0.2 


7.8 ±0.4 


0.07 ±0.02 


0.21 ±0.04 


235.6 ±8.3 


37.7 ±2.5 


-1836.6 ±67.5 


NGC 225 


8.8 ±0.2 


7.4 ±0.2 


0.09 ±0.01 


0.07 ±0.03 


227.7 ±4.5 


33.3 ±0.8 


-1770.9 ± 35.5 


NGC 188 


9.3 ±0.6 


8.1 ±0.8 


0.07 ±0.04 


0.79 ±0.21 


246.7 ± 17.3 


57.1 ±8.8 


-1872.0 ±129.0 


IC 1590 


10.0 ±0.6 


9.2 ±0.6 


0.04 ± 0.02 


0.33 ±0.17 


268.7 ± 17.0 


46.5 ±8.1 


-2092.9 ± 119.8 



Table 3. Age-velocity dispersion relations derived from OCs in 
our present sample with errors in the spatial velocities less than 
20 km s'^ 



Age (Myr) 


o-u 


(TV 


aw 


N 


Age < 500 


28.3 


15.4 


10.6 


339 


500 < Age ^ 1000 


29.1 


13.8 


11.9 


25 


1000 < Age ^ 2000 


31.7 


23.4 


17.7 


16 



4 ORBITAL PARAMETERS AND THEIR 
UNCERTAINTIES 

4.1 The Galactic gravitational potential model 

In this study we employ the axisymmetric Galact ic gravita- 
tional potential model of lAUen fc SantillanI ()l99ll . hereafter 
AS91). This model consists of a spherical central bulge and 
a disc in the form proposed by^Miyamoto & Nagai (1975;), 
plus a massive, spherical halo extending to a radius of 100 
kpc from the centre of the Galaxy. The total mass of the 
model is 9.0 x 10^^ Mq and the local total mass density at 
the solar position is po = 0.15 Mq pc~^. The rotation curve 
of this potential represents the current knowledge of galac- 
tic rotation in the Galaxy. This model is time-independent, 
completely analytical and very simple. The integration of the 
orbit, using this model, is very rapid and can achieve high 
numerical precision. The potential admits two conserved 
quantities, the total energy E and the z-component Jz of 
the angular momentum vector. This model has been used 
to derive the galactic orb i ts of OCs (C94), globular clusters 
llOdenkirchen et al.1 1 19971 : lAllen. Moreno fc Pichaxdol |2006| . 
|2008, hereafter A06, AOS), and clump giants near the Sun 
(SOS). 

4.2 The orbital parameters 

Using the data listed in Table [U the orbits are calculated 
backwards in time over a interval of 5 Gyr. Most clusters in 
our sample have ages less than 100 Myr, they do not even 
move one galactic orbit in the Galaxy. The integration time 
is chosen to ensure clusters can move more galactic orbits 
in the Galaxy and the averaged orbital parameters can be 
determined. For the integration we use the Bulirsch-Stoer 



algorithm of IPress et al.l l|l992l ). The relative change in the 
total energy over the 5 Gyr integration time is of the order 
of 10"^'' to 10-^^ 

The orbital parameters are listed in Tabled Ra and Rp 
are apogalactic and perigalactic distance from the Galac- 
tic centre, which are determined from the averaged max- 
imum and minimum galactocentric distances of the clus- 
ter in the calculated Galactic orbit within the integration 
time of 5 Gyr. The orbital eccentricity e is calculated as 
e = {Ra — Rp)/{Ra+Rp), where Ra and Rp are averages. The 
maximum distance above the Galactic plane, Zmax, is also 
the averaged maximum vertical distances above the Galac- 
tic plane in the cluster's orbit within the given integration 
time. Tp is the orbital period defined as the period of rev- 
olution around the z-axis. Tz is the mean time interval of 
the cluster to cross the Galactic plane from one Zmax to the 
other one in the opposite direction. Because of the right- 
handed orientation of the coordinate system adopted here, 
the negative Jz of the z-component of the angular momen- 
tum vector cor responds to prograde rot ation in the Galaxy 
and vice versa jOdenkirchen et al.lll997l ). 

In Fig.O we plot the distributions of derived orbital pa- 
rameters. In the panels of Zmax, Tp, and Tz, some very large 
values are not plotted. The distributions of Ra, Rp, Tp, and 
Tz are fitted by Gaussian functions ~ g-(^-i^) /2<^ ^ ajj,^ ^jjg 
distributions of e and Zmax are fitted by exponential func- 
tions e~^^^ . The means and dispersions of the Gaussian 
functions for Ra, Rp, Tp, and Tz are derived as: \ir^ — 8.40 
and crj{„ = 1.00 kpc, /ir^ = 7.2S and an^ = 0.96 kpc, 
HTp = 219.3 and ar^ = 22.5 Myr, fir, = 32.3 and ar, = 3.7 
Myr, respectively. The parameters j3 for the distributions of 
e and Zmax are derived as: l3e — O.OS and Pz^a^ = 0.13. The 
derived parameter /3 and Fig.Oshow that, for most of OCs, 
the orbital eccentricities e are less than 0.1 and Zmax less 
than 200 pc. In our sample, the minimum of the Rp is big- 
ger than 1 kpc and the mean of Rp is ~ 7.0 kpc. Therefore, 
the orbits of OCs in our present sample cannot be notice- 
ably affected by a ba rred mass distributi on in the Galactic 
centre within 1 kpc (|Dinescu et all 0*9991 . A06). Fig. [3] also 
indicates that, in one orbital period, most clusters can cross 
the Galactic plane seven times. 

As a representative example, a few of orbits calculated 
with AS91 model in a time-interval of 2 Gyr are presented in 
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Fig. |4l For each cluster, the panel on the left shows the orbit 
projected onto the Galactic plane, while the panel on the 
right shows the meridional orbit. The filled square indicates 
the present observed position for each cluster. The orbits 
of NGC 188, NGC 2682, NGC 2420, NGC 752, and NGC 
2506 are also calculated by C94 and iFinlav et al.l (|l995h . 
For other clusters presented in Fig. |4l Berkeley 33 has the 
maximum of Ra =48.3 kpc, Berkeley 20 has the maximum 
eccentricity e — 0.81, Berkeley 29 and Berkeley 31 have the 
maximum values of Zmax, a nd NGC 67 91 is among the most 
massive OCs known today (|Frielll 19991 ). 

All of the meridional orbits in Fig. 2] are of boxy-like 
type. Clusters move in the meridional plane within the lim- 
ited areas almost filling the boxes symmetrically. But the 
meridional orbits of Berkeley 33 and Berkeley 20 are not 
symmetric with respect to the Galactic plane. The orbits 
projected on the Galactic plane in Fig. |4] indicate the peri- 
odic motions of clusters more clearly. 

4.3 The errors in the orbital parameters 

The errors for the derived orbital parameters listed in Ta- 
blelHinclude two types of uncertainties affecting the derived 
results. The first one is the intrinsic variation of the orbital 
parameters within the 5 Gyr integration interval. The dis- 
persions of the averaged orbital parameters are calculated 
over the number of galactic orbits, which indicate the in- 
trinsic nature of the orbit, and this may be due to effects 
from chaos, and/or a complex distribution of orbit families 
l|Dinescu et alll 19991 ). 

On the other hand, the main errors of the orbital param- 
eters come from the observational uncertainties of the input 
data. In the input data, a 20% relative error is assumed for 
distance, the median of the relative errors in radial veloc- 
ities is about 6% and the median of the relative errors in 
proper motions is 23%. The uncertainties in distance and 
proper motions are the main sources for the errors in the 
derived orbital parameters. The effects of observational un- 
certainties cannot be simply propagated into the derived or- 
bital paramet ers |Od cnkirchcn fc Brosche 1992 ). Following 
iDinescu et all (|l999l ). the initial conditions are generated in 
a Monte Carlo fashion by adding Gaussian deviates to the 
observed absolute proper motions, radial velocities, and dis- 
tances. The standard deviations are taken to be the errors 
for the input data listed in Table [T] For each cluster, the 
errors for its orbital parameters have been calculated based 
on 1000 separate integrations. 

Compared with the uncertainties in the derived orbital 
parameters due to the errors in the input data, the intrinsic 
uncertainties of the orbits within the given integration in- 
terval are very small. For each cluster, only considering the 
observational errors of the input data, the relative errors in 
the derived orbital parameters are also calculated. For each 
orbital parameter, the median of the relative errors for this 
parameter is calculated and listed as percentage in column 
2 of Table [5] Table [5] shows that the uncertainties of orbital 
eccentricity e and the maximum distance from the Galac- 
tic plane ^max are much more affected by the observational 
errors. The large uncertainty of e is propagated from the 
uncertainties of both Ra and Rp. The reason for the large 
uncertainty of Zmax is that, for most OCs, they move near 
the Galactic disc and have small vertical distances from the 



Table 5. The relative error percentages of derived orbital param- 
eters due to observational uncertainties and the relative difference 
percentages of derived orbital parameters by different Galactic 
potential models. 



Parameter observational uncertainty P90 FSC96 

(per cent) (per cent) (per cent) 



Ra 2.9 0.0 0.2 

Rp 3.5 0.5 0.3 

e 27.5 0.0 0.0 

2max 43.3 1.1 21.2 

Tp 2.8 0.8 1.5 

3.9 2.9 31.5 



Galactic plane, so small changes in input data produce big 
relative changes of their orbits in the direction perpendicular 
to the disc. 

4.4 The effect of adopting different observational 
errors in input data on the derived orbital 
parameters 

Comparing the observational errors in the input data, the 
relative error in radial velocity is smaller than that in proper 
motion and distance. The errors in the derived orbital pa- 
rameters are dominated by the observational errors in proper 
motions and distances. In order to estimate the effect of 
different observational accuracy in proper motions and dis- 
tances on the derived orbital parameters, we repeated the 
Monte Carlo procedure mentioned above. We assumed 50% 
and 80% relative errors for the adopted proper motion and 
distance data, respectively. We also assumed 50% and 80% 
relative errors for both proper motion and distance all to- 
gether. 

There are six sets of simulations. In each set of simu- 
lation, except for the assumed different relative errors for 
proper motions or distances, the errors for other input data 
are the same as those considered in the previous section. For 
each set of simulations, we calculated the relative error for 
each orbital parameter, the medians of the relative errors 
are listed in Table [6] Table |6] indicates that the relative er- 
rors in e and Zmax increase obviously with the increase of 
errors in input data. The uncertainties of e and Zmax are 
very sensitive to the errors of the input data just as Table [5] 
has indicated. Table [S] also shows that, assuming the same 
relative error for proper motion and distance, the relative 
errors in the derived orbital parameters caused by the er- 
ror of distance are bigger than those caused by the error of 
proper motion. 



5 ORBITAL PARAMETERS IN DIFFERENT 
GALACTIC POTENTIAL MODELS 

We use two different potential models to calculate the dif- 
ferences in the derived orbital parameters. The first model 
is the one pro posed by iPacz yrlski (1990, hereafter P90) and 
then used by iDinescu et all |l999j) as their representative 
model to calculate the orbits of 38 globular clusters. This 
model consists of axisymmetric potential with three compo- 
nents: bulge, disc, and dark halo. The bulge is modeled as 
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Table 6. The relative error percentages of derived orbital parameters for adopted different observational uncertainties in proper motions 
fj, and distances from the Sun dQ . 



Parameter observational uncertainty 





50% in ^l 


80% in n 


50% in dQ 


80% in dQ 


50% in /x+d0 


80% in ^l+dQ 




3.9 


4.8 


6.8 


11.5 


7.5 


14.5 


Rp 


5.0 


5.8 


6.7 


9.9 


7.7 


12.6 


e 


45.0 


54.4 


45.0 


62.5 


61.4 


102.2 




76.8 


112.9 


81.5 


130.0 


121.0 


271.7 


Tp 


3.7 


4.3 


5.8 


9.7 


6.5 


12.3 




6.3 


11.7 


8.3 


14.9 


11.1 


26.2 



Plummer potential (IPlummei 191ll). The disc is the same 
as AS91 in the form of lMivamoto fc Nagail l| 19751 ) with dif- 
ferent coefficients. The dark halo is modeled as logarithmic 
potential, which assures a flat rotation curve, but imply an 
infinite mass. 

The second model is the on e pr oposed by 
iFlvnn. Sommer-Larsen fc ChristensenI (|l996l . hereafter 
FSC96) and used by N04 to calculate the orbits of nearby 
F and G dwarf stars. This axisymmetric potential also 
consists of three components: central core, disc, and dark 
halo. The potential of the central core is modeled by two 
spherical components, representing the bulge/stellar-halo 
and an inner core component. The disc is modeled using a 
combi nation of three analytical discs of lMivamoto fc Nagail 
The potential of the dark halo is assumed to be 
logarithmic. 

In Figs. [S] and [S] we compare the orbital parameters de- 
rived from P90 and FSC96 models with those derived from 
AS91 model. Fig. [5] shows that the orbital parameters de- 
rived from P90 model are very consistent with those de- 
rived from AS91 model, and no systematic differences can 
be found. For FSC96 model. Fig. |6] shows that Rp and e 
are consistent between FSC96 and AS91 models, and there 
are no systematic differences in these two parameters. But 
for Ra > 15 kpc, the derived Ra and Tp for FSC96 model 
are systematically smaller than those derived from AS91 
model. Fig. |5] also indicates that the derived 2:max and 
for FSC96 model are systematically bigger than those de- 
rived from AS91 model. 

In Figs. [7] and [H for P90 and FSC96 models, we show 
the orbits for the same clusters as presented in Fig. 3] Fig. [7] 
indicates that the orbital shapes for most of clusters, in- 
cluding the asymmetric meridional orbits of Berkeley 33 and 
Berkeley 20, are very similar between P90 and AS91 mod- 
els. Fig.[8]shows that for most of clusters, the orbital shapes 
derived from FSC96 model are different from those derived 
from AS91 model. The meridional orbits of Berkeley 33 and 
Berkeley 20 are symmetric in FSC96 model. 

For each cluster, the relative difference in each or- 
bital parameter from different models is calculated as 
(-P(P90 or FSC96) " 'PaS9i)/-PaS91' ^^ere Ppgg, ^'aS91' 
and -Ppscgg represent the derived orbital parameters by 
P90, AS91, and FSC96 models respectively. For each orbital 
parameter, the median of the relative difference between the 
given models and AS91 model is listed in Table [5] Table [5] 
shows that the relative differences of the derived orbital pa- 
rameters for P90 model are smaller than those for FSC96 
model. Especially for Zmax and Tz, the results derived from 
FSC96 model are very different from those derived from 



AS91 and P90 models. The consistent orbital parameters 
derived from AS91 and P90 models are due to the very sim- 
ilar mass distributions of these two models. The main dif- 
ference between AS91 and P90 models is the bulge model, 
whereas bulge model dominates the orbit of a cluster only 
within 1 kpc from the Galactic centre. No clusters in our 
present sample can move to this range. The FSC96 model 
is very different from AS91, the total mass interior to the 
same galactocentric distance R in the FSC96 model is much 
larger than that in the AS91 model when R > 10 kpc. The 
larger mass of FSC96 model can control the cluster to move 
within smaller apogalactic distance. On the other hand, the 
mass of disc in FSC96 model is smaller than that in AS91 
model, clusters can move to more larger vertical distances 
in FSC96 model. 

Table [5] also shows that, for most of derived orbital pa- 
rameters, the relative differences due to different potential 
models are smaller than those from observational errors. But 
for FSC96 model, the relative differences in Tz are much 
larger than the uncertainties due to observational errors of 
input data, which indicates the major effect of the disc model 
of FSC96. 



CHARACTERISTICS OF THE ORBITAL 
PARAMETERS 



6.1 Relations between Ra 

and Zmax — \z\. 



The left panels of Fig. [9] show the Ra vs Rqq, Rp vs Rqq, 
and Zmax vs |z| diagrams, where 1^1 is the current observed 
vertical distance from the Galactic disc. The panel of Ra vs 
Rqq diagram indicates that, for most of clusters in present 
sample, their current observed positions are very close to 
their apogalacticons. Whereas, the panel of Rp vs Rqq dia- 
gram indicates that most of clusters are far away from their 
perigalacticons. The panel of 2max vs diagrams indicates 
that most of clusters do not arrive at their maximum dis- 
tances from the disc and are crossing the Galactic plane. 

In the right panels of Fig. (9] we plot the histograms of 
the relative differences for corresponding parameters in the 
left panels. The top panel shows the relative difference be- 
tween Ra and Rqq and indicates that, at present, about 
70% clusters are moving within only about 5% distances 
from their apogalacticons. The median of the relative dif- 
ferences between Ra and Rqq is 2.0%. The middle panel 
shows the relative difference between Rp and Rqq, and in- 
dicates a significantly larger relative difference with respect 
to that between Ra and Rqq- The median of the relative 
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differences between Rp and Rqq is 8.1% which is larger 
than that between Ra and Rqq- The bottom panel shows 
the relative differences between Zmax and \z\, and indicates 
very large differences between these two parameters. The 
median of the relative differences is 68.8%. The mean time 
Tz for our present sample is very short, so the time for a 
cluster moving near its Zmax is short, which is the possible 
reason for the large differences between the observed z and 
the orbital parameter Zmax- 

For two clusters NGC 2682 and NGC 2420 whose or- 
bital eccentricities are larger than 0.1, C94 calculated the 
probability of finding the clusters at the given galactocen- 
tric distances and found that the detection probability is the 
largest at the cluster's apogalacticon. Following C94, we also 
define a probability function P{R) oc -^j^y to calculate the 
probability of a cluster to be observed at the galactocentric 
position R during the orbital period Tp, where ^{R) is the 
cluster's velocity at R. Fig. 1101 shows the probability distri- 
butions P{R) for the same clusters as in Fig. |4l The filled 
square in Fig. [TO] indicates the present observed position 
for each cluster. Fig. [10] shows that the detection probabil- 
ity P{R) increases with the galactocentric distance R, the 
largest detection probability is at the cluster's apogalacti- 
con. The derived probabilities for other clusters also indicate 
that the detection probability for a cluster at its apogalac- 
ticon is the largest one during its orbital period, which is 
consistent with the result of C94. Because the detection 
probability at the apogalacticon for a cluster is the largest 
during its orbital period, it is easier to find a cluster near 
its apogalacticon just as Fig. [9] indicates. 



6.2 The radial metallicity gradient 

The OCs have long been used as tracers of radial metallic- 
ity gr a dients in the Galactic disc. Since the early work by 
IJanesI (|l979l ). others have found general agreement in the 
existence and magnitude of the trend. Most investigators 
have found gradients of —0.06 to —0.09 dex kpc~^ over a 
range of distances from 7 to 16 kpc from the Galactic c entre 
iTwarog. Ashma fc Anthonv-Twarodl 19971 : iFYielll 19991 ). Re- 



cently, Friel et al] 1 20021 ) presented metallicities for a sam- 
ple of 39 intermediate age and old OCs based on an up- 
dated abundance calibration of spectroscopic indices. They 
found a metallicity gradient of — 0.06±0.01 dex kpc~^ over a 
range of galactocentric distances of 7 to 16 kpc. I Chen et al.l 
(|2003l ) compiled a OCs catalogue of 119 objects with ages, 
distances, and metallicities available, which led to a metallic- 
ity gradien t of — 0.063± 0.008 dex kpc~^, similar to the result 
derived bv lFriel et al. I (2002i ) from a homogeneous sample. 

The [Fe/H] vs current observed galactocentric distances 
Rqq diagram is plotted in the top panel of Fig. 1111 109 
clusters with [Fe/H] data listed in our present sample are 
plotted as open circles, 48 clusters not listed in our present 
sample but with [Fe/H] data are plotted as plus signs. There 
are 12 clusters in the range of Rqq > 13.5 kpc, only 3 
clusters are included in our present sample. Considering the 
observational errors in [Fe/H] and Rqq, we perform a linear 
least-square fitting to the clusters listed in our sample with 
Rqq < 13.5 kpc. For each cluster, we assign the typical 
observational uncertainty of 0.15 dex to their [Fe/H] data. 
The errors of Rqq are calculated from the data listed in 



Table [2] We get a gradient of -0.070 ± 0.011 dex kpc"\ 
We also perform the same fitting to all clusters in the range 
of Rqq < 13.5 kpc, we get a gradient of -0.069 ± 0.008 
dex kpc~^. The fitted straight-line for the clusters listed in 
our present sample within Rqq < 13.5 kpc is plotted in 
the same panel of Fig. [TT] and the slope b of this line is also 
labeled. 

In the most recent study, IChen et al.l l|2007l ) derived 
a radial metallicity gradient of —0.058 ± 0.006 dex kpc~^ 
based on a sample of 144 OCs. Our derived radial metal- 
licity gradi ent of —0.0 7 is slightly smaller than the previ- 
ous results (|Friel et al.| [2002: Chen et al. 2003, 20"07l). I3ased 
on 45 OCs with high-resolution spectroscopv. iMagrini et al.l 
(|2009l ) found a steep metallicity gradient for clusters with 
Rqq < 12.0 kpc and a plateau for clusters at larger galac- 
tocentric distances. Fig.[TT]indicates that the metallicity dis- 
tribution of clusters with Rqq > 13.5 kpc is flat and there is 
no signiflcant radial metallicity gradient withi n this distance 
range, which is consistent with the result of lMagrini et al.l 
(2009). In our present sample, only three OCs have galacto- 
centric distances greater than 13.5 kpc (the open circles in 
Fig. Ilip . In order to diminish the effect of the small sample 
at large galactocentric distances, we only consider OCs with 
Rqq < 13.5 kpc in our present sample to derive the radial 
metallicity gradient. The maximum galactocentric distances 
of the OCs in the previous studies extend to 17 kpc and some 
clusters within t he fiat metallicity distribution were include d 
in their samples dFriel et al.l |200^ : IChen et all |2003| . |2007| ). 
which is the main reason why they obtained larger radial 
metallicity gradient. If we use all clusters with Rnn < 17.0 
kpc, we get a gradient of —0.056 ± 0.007 dex kpc~ which i s 
consistent with the most recent result of lChen et al.l |20o3) ■ 

The radial metallicity gradient of the OCs provides 
strong constraints on the formation and evolution of the 
Galaxy. Detailed models of Galactic chemical evolution have 
been improved over the last decades and most models can 
reproduce the presentl y observed radial m etallicity distri- 
bution of the Galaxv dMagrini et al. l2009l . and references 
therein). More recently. IFu et al.l 12009') considered various 
mechanisms including infall, star formation and delayed disk 
formation to find the effect of each mechanism on their de- 
rived galactic chemical evolution model. They found that 
using the star formation rate(SFR) of the modified Kenni- 
cutt law, their model can properly predict both the current 
metallicity gradient and its time evolution. But, their best 
model also predicts that the outer disc has a steeper gradient 
than the i nner disc, whic h is co ntrary to the result derived 
from OCs. iMagrini et al. I l|2009l ) adopted an inside-out for- 
mation model of the Galactic disc to reproduce the radial 
metallicity gradient of the OCs. In their model, the infall 
of gas is represented by an exponential law combined with 
the distribution of gas in the halo. The inner regions are 
rapidly evolving due to the higher infall and SFR, while the 
outer parts evolve more slowly. Their model can reproduce 
the main features of the metallicity gradient and the evolu- 
tion of the OCs. In order to better reproduce the metallicity 
plateau at large galactocentric distances, an additional uni- 
form infiow per unit disc area should be considered, but it is 
difficult to reconcile with the present-day radial distribution 
of the SFR. A sequence of merging episodes in the past 
history of the Gal axy would explain t he outer metallicity 
plateau of the OCs (|Magrini et al.ll2009l ). 
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In our previous discussion, we have showed that the 
apogalacticon is the place where a cluster spends the largest 
fraction of its life. In the bottom panel of Fig. 1111 we plot the 
[Fe/H] vs Ra diagram for clusters listed in our sample within 
Rqq, < 13.5 kpc. We also perform a linear least-square 
fitting to the data considering the observational errors in 
[Fe/H] and Ra- The errors of Ra are taken from Table |3l 
We plot the fitted straight-line in the same panel and also 
label the slope h = —0.082 ± 0.014 dex kpc~^ in the panel. 
This result indicates that, for clusters within Rqq < 13.5 
kpc, the observed metallicity gradient at present is similar 
to that derived from the most probable observed positions 
of the clusters, which is consistent with the result of C94. 



6.3 Comparison of orbital eccentricities for 
different populations 

Fig. [12] shows the histograms of orbital eccentricities for dif- 
ferent populations: globular clusters (top left panel), disc 
giants (top right panel) , disc F and G dwarf stars (lower left 
panel), and OCs (lower right panel). Each histogram is 
normalized to its maximum value in the distribution. The 
orbital parameters for 54 globular clusters calculated with 
AS91 model are taken from A06 and AOS. The orbital pa- 
rameters of disc clump giants calculated with AS91 model 
are taken from SOS. The orbital parameters of disc F and 
G dwarf stars calculated with FSC96 model are taken from 
N04. 

Globular clusters represent the halo population which 
is primarily a system with a large velocity dis persion, and a 
wide range of orbit characteristics is expected (jPinescu et all 
This is indeed what the top left panel of Fig. ll2l shows. 
There is a large range in orbital eccentricities of globular 
clusters, which can get as low as 0.1, but, in the mean, the 
eccentricities are high with an average of ~ 0.5. 

SOS assigns to each star its probability to belong to the 
thin disc, the thick disc, the Hercules stream and the halo 
based on the basis of its ([/, V, W) velocity and their veloc- 
ity ellipsoids. In the top right panel of Fig. 1121 we show the 
histograms of orbital eccentricities for stars with probabili- 
ties higher than S0% to belong to the thin disc (solid lines) 
and stars with probabilities higher than 80% to belong to 
the thick disc (dot lines). Fig. [T2] indicates that no thin disc 
giant has eccentricity greater than 0.3 and no thick disc gi- 
ant has eccentricity less than 0.1. The orbital eccentricities 
of the thin and the thick disc giants are overlapped between 
0.1 and 0.3. The average of orbital eccentricities for the thin 
disc giants is ~ 0.1, and the one for the thick disc giants is 
~ 0.4, which is less than the average of orbital eccentricities 
for globular clusters. 

The histogram of orbital eccentricities for disc F and 
G dwarf stars in the lower left panel of Fig. [12] shows that 
most of stars belong to the thin disc whereas a fraction of 
stars belong to the thick disc. If we assume that stars with 
eccentricities greater than 0.3 belong to the thick disc, the 
thick-disc fraction i s about 3.6% wh i ch is c lose to the value 
of 2.9% derived by iHolmberg et all (120071 ) from the same 
sample with different method. 

The lower right panel of Fig. [12] shows the histograms 
of orbital eccentricities for OCs in our present sample calcu- 
lated with AS91 model (solid lines) and FSC96 model (dot 
lines). The distributions of eccentricities in these two mod- 



els are very similar. The distribution of orbital eccentricities 
for OCs is similar to that of disc F and G dwarf stars which 
indicates that most of clusters belong to the thin disc and 
a fraction of thick-disc clusters exist in our sample. We as- 
sume the same limited value of 0.3 for orbital eccentricities 
to distinguish thin-disc and thick-disc clusters, and find 3.7% 
clusters in our sample are probably thick-disc clusters. 

In our present sample, the Zmax values of Berkeley 29, 
Berkeley 31, and Berkeley 33 are 15.1 ± 15.4, 16.2 ± 12.2, 
and 8.5 ± 9.3 kpc, which are much larger than those of the 
other clusters in the sample. The maximum of Zmax of the 
other clusters is less than 3 kpc and the observed maximum 
of z is 2.1 kpc for Berkeley 29. The very large errors of Zmax 
for these three clusters indicate that their orbits are very 
uncertain and more observations are needed to improve the 
precision of the input data for these clusters. Berkeley 20 
has the maximum of e = 0.81 ± 0.17 in our present sample 
which indicates it is a halo cluster. But this cluster has an 
age of 6.0 Gyr and [Fe/H]= -0.61 ± 0.14 indicating it is 
a thick disc cluster. If we adopt the limited value e — 0.3 
with 3a from the mean, the cluster can be identified as a 
member of the thick disc population only based on orbital 
eccentricity e. More observations are also needed to improve 
the orbital parameter for this cluster. 

6.4 Comparison with previous results and the 
open cluster complexes 

There are two studies which calculated the orbits for a large 
samples of OCs and the orbital parameters were also de- 
rived. BKO calculated the orbits for 69 OCs and the orbital 
parameters Ra, Rp, e, and 2max were also derived for those 
clusters. In Fig. 1131 for common clusters, we compare the 
derived orbital parameters in AS91 model with those de- 
rived by BKO. Fig. 1 131 indicates that, for most clusters, the 
derived Ra, Rp, and e are consistent with those derived by 
KBO. The maximum difference in Ra is for NGC 2420. The 
radial velocity for this cluster is 73.6 km s~^ in our present 
study, but the value of 115 km s~^ was adopted by BKO. 
The large difference in radial velocity make the large differ- 
ence in Ra for NGC 2420. The maximum difference in Rp 
and e is for NGC 7789. The radial velocity of -64 km s"^ 
is adopted for this cluster in our study, BKO adopted the 
value of —32 km s~^. The large difference in Rp and e for 
NGC 7789 is also due to the large difference in radial veloc- 
ity data. For most of OCs, the orbital parameter Zmax in our 
present study are less than those derived by BKO. The sys- 
tematic difference in Zmax is due to the different disc models 
ado pted by this stu d y and by BKO. 

jPiskunov et al.l ( 20061 ) calculated the orbits for a sam- 
ple of 148 OCs within d^y ^ 0.85 kpc and the mean or- 
bital parameters are listed in their Table 2: — 8.631 
kpc, — 6.706 kpc, fi^ = 0.127, and — 0.260 kpc. 

For the clusters in our present sample within the same dis- 
tance range, the corresponding mean orbital parameters cal- 
culated with AS91 model are derived as: /i_r^ — 8.289 kpc, 
fiRp = 7.378 kpc, fie = 0.059, and = 0.084 kpc. Be- 

cause the input data for most of clusters in this stud y are 
similar to those in the study of jPiskunov et al.l (|2006l ). the 
differences of the derived orbital parameters are mainly due 
to the adopted different Galacti c models. jPiskunov et al] 
(j2006l ) used the Galactic model of lSaio fc Yoshiil l| 19791 ). the 
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analytic solution of the orbital parameter can be derived 
from their model. In comparison, the much larger disc mass 
in AS91 model is responsible for the large difference in the 
derived value of Zmax with respect to the other model. 

Many surveys have found that the structure of the 
Galaxy is more complex than previously thought, and about 
10 moving groups can be identified from nearby stars wit h 
heliocentric distances less than 100 pc ()Bovv et al.ll2009l ). 
Using OCs with the projected distance onto the Galactic 
plane d^y ^ 0.85 kpc from the Sun, and based on the surface 
density distributi on and th e tangential velocity distribution 
of those clusters, Pisk unov et al. (2006) identified four open 
cluster complexes (OCCs). The numbers of kinematic mem- 
ber clusters for the four OCC s: OCC 1, OCC 2, OC C 3, and 
OCC 4 are 23, 27, 8, and 9 ijPiskunov et al.ll2006l l. In our 
present sample, for the corresponding OCCs, t he numbers of 
member clusters identified by Pisku nov et al.l (|2006) whose 
orbits can be determined are 20, 13, 3, and 4. 

Those OCCs were detected from their overdensity in the 
spatial distribution of OCs in the Solar neighbourhood and 
the membership of the member clusters were determined 
only based on their tangential velocities ()Piskunov et all 
I2OO6I ). In Fig. [H we use the Vac — Rac and z^ax — e 
diagrams to attempt to recover those OCCs, where Vac 
is the total galactocentric velocity of cluster. The crosses 
in Fig. [14] represent clusters with heliocentric distance less 
than 1.3 kpc, the open circles represent member clusters in 
each OCC. In the panels of OCC 1, the open circles rep- 
resent member clusters with age less than 30 Myr and tri- 
angles represent member clusters with age between 30 and 
80 Myr. There are 18 member clusters in OCC 1 whose age 
data are available. The Vgc — Rgc diagram of OOC 1 in- 
dicates that 6 clusters with age less 30 Myr and 4 clusters 
with age great than 30 Myr distribute in a small region with 
212 < Vac < 219 km s~^ and 8.0 < Rac < 8.5 kpc, those 
clusters can be considered as kinematic members of OCC 1. 
But in the Zmax — e diagram of OCC 1, it is difficult to 
identify any clustering of those member clusters. We could 
only determine the orbits of 50% of the assumed 
members of OCC 2. Most of member clusters of OCC 2 
have age between 200 and 400 Myr. The panels of OCC 2 
indicate that we cannot detect any clustering of the member 
clusters in the Vac — Rgc and Zmax — e diagrams. The 
OCC 3 and 4 only include few of member clusters, the pan- 
els of those two OCCs indicate no clustering of the member 
clusters can be identified. Fig. [14] indicates that only young 
OCCs can be identified in the Vgc — Rcc diagram, and it 
is difficult to detect clustering in OCs with the Zmax — e 
diagram. 

More recently, using the radial velocity, proper mo- 
tion, inclination, and Galactic latitude of 341 OCs 
with age les s than 100 Myr and within 2.5 kpc 
from the Sun, Ide la Fuente Marcos fc de la Fuente Marcos! 
(|2008t ) studied the clustering in those clusters. Most 
of member clusters in the closest OCC detected by 
Ide la Fuente Marcos fc de la Fuente Marcos! (|2008! ') have 
age less than 30 Myr, which is consistent with that we 
have found in the Vgc — Rcc diagra r n of OCC 1. 
!de la Fuente Marcos fc de la Fuente Marcos! (I2OO8! ') pointed 
out that 20 Myr is the characteristic timescale of an OCC, 
after this timescale, the complex may no longer be recogniz- 
able in the space of the orbital elements as the majority of its 



member have evaporated; only in the context of corotation 
resonances or mergers, the nongenetically related dynamical 
groups of old OCs can be detected in the Galactic disc. 



7 SUMMARY AND CONCLUSIONS 

We have presented a sample of 488 OCs whose distances, 
radial velocities, and absolute proper motions are used to 
derive their orbits in the Galaxy. The kinematical and the 
orbital characteristics of this sample are analyzed. The main 
results are listed as following: 

• For OCs with errors in the spatial velocities less than 
20 km s~^ in our present sample, the velocity ellipsoids are 
derived as (ay, ay, o-w ) =(28.7, 15.8, 11.0) km s The 
ages for most of clusters in our present sample are less than 
500 Myr, and this sample represents a young thin disc popu- 
lation. The velocity dispersions of OCs in the three velocity 
components increase with the age of the cluster subsample, 
which indicates the continuous heating of the disc (N04). 

• The orbits of OCs are calculated with three Galac- 
tic gravitational potential models. Considering the intrinsic 
variation of orbital parameters and the effects of observa- 
tional uncertainties, the errors for the derived orbital pa- 
rameters are determined with a Monte Carlo fashion. The 
major errors come from the observational uncertainties of 
the input data and are mainly affected by the errors in the 
distance data. The observational uncertainties mainly affect 
the derived orbital eccentricities e and the maximum dis- 
tances above the Galactic disc Zmax. 

• The total masses within given distance range for differ- 
ent Galactic models can affect the derived orbital parameter 
Ra and the orbital period Tp. The disc models for different 
Galactic models mainly affect the vertical movement of clus- 
ters and change the derived Zmax and T^. The uncertainties 
in the derived orbital parameters caused by different models 
are smaller than those caused by observational errors in the 
input data. 

• The detection probability for a cluster at the given 
galactocentric distance is calculated and the largest detec- 
tion probability is at the cluster's apogalacticon. For most of 
OCs in our present sample, their present observed positions 
are very close to their apogalacticons. The mean of orbital 
period for OCs is about seven times longer than the time 
for clusters crossing the Galactic plane. 

• Based on the presently observed galactocentric dis- 
tances, the radial metallicity gradient for clusters with 
Rqq < 13.5 kpc is derived with a slope —0.070 ± 0.011 
dex kpc~^, whic h is consistent wi th the previous studies 
(|Friel et al.!!2002! : !Chen et al.l!200j ). The radial metallicity 
gradient derived based on apogalactic distances for the same 
sample is -0.082 ± 0.014 dex kpc"^ 

• The orbital eccentricities e for different populations: 
globular clusters, disc giants, disc F and G dwarf stars, and 
OCs are compared. The orbital eccentricities for globular 
clusters occupy a large range and with a mean of ~ 0.5. 
The mean of orbital eccentricities for the thick disc giants 
is ~ 0.4 which is bigger than the mean of ~ 0.1 for the thin 
disc giants, F and G dwarfs, and OCs. There are about 3.7% 
clusters in our sample belonging to the thick disc. 

• Using the Vgc — Rgc diagram, only one OCC 
could be identified and most of its members are 
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younger than 30 Myr. We find it difficult to detect 
any OCCs in the Zmax — e diagram. 
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Figure 1. The volume density distributions of OCs as a function of their heUocentric distance cIq. The open circles represent the OCs 
in DAML catalogue with distance data available, the filled circles represent the OCs in our present sample. 
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ure 2. Distributions of ages, metallicities [Fe/H], and observed galactocentric distances Rqq of OCs in our present sample. 
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Figure 3. Distributions of derived orbital parameters calculated with AS91 model for OCs in our present sample. 
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Figure 4. Meridional Galactic orbits and orbits projected onto the Galactic plane in the time-interval of 2 Gyr for some OCs calculated 
with AS91 model. The filled square shows the present observed position for each cluster. 
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Figure 5. Comparison of derived orbital parameters with AS91 and P90 model. 
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Figure 6. Comparison of derived orbital parameters with AS91 and FSC96 model. 
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Figure 7. Meridional Galactic orbits and orbits projected onto the Galactic plane in the time-interval of 2 Gyr for some DCs calculated 
with P90 model. The filled square shows the present observed position for each cluster. 
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Figure 8. Meridional Galactic orbits and orbits projected onto the Galactic jjlaiic in the time-interval of five times orbital periods 
some OCs calculated with FSC96 model. The filled square shows the present observed position for each cluster. 
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Figure 9. Relations between the orbital parameters Ra, Rp, ^max of OCs in our present sample and their observed data Rqc, \z\- 
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Figure 10. The detection probability distributions P{R) for some OCs in our present sample.The filled square is the present observed 
position for each cluster. 
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Figure 11. Metallicity gradients for OCs. Clusters included in our present sample are plotted as open circles, clusters not included in 
our present sample but with [Fe/H] data are plotted as plus signs. Top panel: The metallicity gradient of OCs is derived based on the 
currently observed galactocentric distances of OCs Rqq ; bottom panel: The metallicity gradient of OCs is derived based on apogalactic 
distances Ra- The straight-lines are the best-fitting results of the metallicity gradients for OCs in our present sample with Rqq < 13.5 
kpc and the slopes b for each line, i.e. the gradients, are labeled in each panel. 
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Figure 12. Comparison of orbital eccentricities e for different populations: globular clusters (A06+A08), disc giants (SOS), disc F and 

G dwarf stars (N04), and OCs (this work). 
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Figure 13. The derived orbital pajrameters Ra, Rp, e, and Zmax with AS91 model in this work are compared with those derived by 
BKO. 
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Figure 14. The orbital parameter distributions of four open cluster complexes (OCCs) identified bv lPiskunov et al] | |200(J) . Open circles 
are member clusters for each OCC and crosses are clusters with heliocentric distance less than 1.3 kpc. In the panels of OCC 1, open 
circles represent member clusters younger than 30 Myr and triangles represent member clusters with age between 30 and 80 Myr. 



